Two-photon width of r\ c and rf c from Heavy-Quark Spin Symmetry 
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We evaluate the two-photon width of the pseudoscalar charmonia, T~ ll {ric) and T 11 {ri' c ), within a 
Heavy-Quark Spin- Symmetry setting and show that whereas the former width agrees with experi- 
ment, the latter is more than twice larger than the recent measurement by CLEO. When binding- 
energy effects are included in the rj c case, the discrepancy is worse, pointing out at a possible 
anomaly in the rj' c decay. 
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I. INTRODUCTION 

Whereas heavy-quarkonium production is still a great 
source of debates (see 0, for recent reviews), the 
physics of quarkonium decay seems to be better under- 
stood within the conventional framework of QCD. How- 
ever, a recent estimation of the ratio of the two-photon 
width of the rf c to that of the r\ c by the CLEO collabora- 
tion Q seems to contradict most of the theoretical pre- 
dictions HM. Indeed, by assuming B(r\ c — > KKtt) = 
B(rf c — > KKtt), they have obtained T^Jrj' c ) = 1.3 ± 0.6 
keV, whereas the predictions of 0, 0, |(J range from 3.7 
to 5.7 keV. 

It is our purpose here to have another look at this 
problem using an effective Lagrangian procedure satis- 
fying heavy-quark spin symmetry and including binding 
energy or equally mass effects, which would take into ac- 
count features typical of radially-excited states. 

Indeed, the r]' c is the first radially-excited pseudoscalar 
charmonium, labeled in the spectroscopic notation by 
2 1 S , with a mass M„/ = 3638±5 MeV Q, that is notice- 
ably higher than for the ground states rj c and J/ip. Its 
first observation was done by the Belle collaboration j|j 
in B — > KKsK~tt + decay and was further confirmed by 
BaBar 0| . 

As a consequence of C-conservation, rj' c , like rj c , can de- 
cay into two photons, which is from a theoretical point of 
view a rather clean channel to analyze. There have been 
several calculations of the rf c — > 77 in the literature, some 
following Bethe-Salpeter equation [Io|. follo wing Salpeter 
equation or relativistic quark models 0, yj [ll|, [lj , and 
some based on the non-relativistic results (see for in- 
stance 0|) but taking into account differences in the 
singlet and triplet wave function at the origin 0, [w| . In 
particular, non-relativistic calculations can only be done 
by considering the rj' c as a 2 S state with the same mass, 
2m c , as the IS state with the result that the calculated 
decay rate differs from the r\ c one only through the wave 
function at the origin, see Eq (3.17) of [13|, or the long 
distance matrix element of NRQCD, see Eqs. (4.17) and 
(4.19) of 0. 



Since the i)' c is more than 600 MeV above the 77 c , the 
mass effects on the decay rate could be important. A 
better approach, which would allow the inclusion of such 
effects, would be to use relativistic kinematics in the cal- 
culation of the width. For this purpose, we need to con- 
struct an effective Lagrangian for the process cc — > 77 
by expanding the charm-quark propagator in powers of 
q 2 /m 2 , with q = p c — p s , and neglecting terms of 
0(q 2 /m 2 ) terms. The propagator will now depend only 
on the charm-quark mass and the binding energy of the 
charmonium state 0, 0] . 

The effective Lagrangian derived in our approach will 
then allow a calculation of the decay amplitude in terms 
of the matrix element of a local operator. The latter is, 
for the two-photon decay width of rj c and rf c , the matrix 
element for the axial-vector current 07^750 between the 
vacuum and ij c or rj' c . 

The non-perturbative parameters are here the decay 
constant /„i s and f n 3 g 1 which can be given by the spa- 
tial wave function at the origin ip(0) |20|. It should be 
stressed here that our approach differs from the tradi- 
tional approach in an important way. We express the 
decay amplitude in terms of the matrix element of a lo- 
cal operator which could be measured or extracted from 
measured physical quantities, like the leptonic-decay con- 
stant or could also be computed via sum rules |2l| or 
lattice simulations [22^. 

We shall rely on the heavy-quark spin-symmetry 
(HQSS) relations 0, 0, which state the equality 
between f Va and /ju and between f v > and fy/. The 
derivation of these relations is based on the fact that, 
in our approach, the flavor-conserving charm-quark cur- 
rents c"7 M c and 07^,75 c take the form of an effective cur- 
rent in which c and c are replaced by static heavy quark 
field operator and the C(l/(2m c )) terms are neglected. 

In this paper, we shall first derive the effective La- 
grangian for the decay of singlet IS state of charmonium 
into two photons. We then show that this effective La- 
grangian, combined with HQSS, gives the same result as 
the traditional non-relativistic approach and produces a 
decay rate for rj c in agreement with measurement. In the 
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next section, we use this Lagrangian to determine the rf c 
two-photon decay rate in terms of the ip' leptonic width, 
as our main purpose here is to see whether our approach, 
which works for the rj c , could explain the observed decay 
rate of the rj' c into two photons. 



II. EFFECTIVE LAGRANGIAN FOR 1 S 
DECAY INTO TWO PHOTONS 

As announced, we now write down an effective La- 
grangian for the coupling of the cc pair to two photons 
and to a dilepton pair 11: 



(1) 



with c\ 



and C2 = 



c(47ra em ) 



Ml + b^M^ 
The factor lJJM^+b Va M T)< .) in c\ contains the binding- 



energy effect [1$, |ljj an d is obtained from the denomi- 
nator of the charm-quark propagator (fci , k 2 being the 
outgoing-photon momenta): 



1 



[(kx - k 2 ) 2 /i - m 2 



(2) 



by neglecting the term containing the relative momenta 
q = p c — p s of the quarks. For real photons, this factor 
can be written as 



1 



[(M 2 + bM)/2] 



(3) 



with b (= 2m c — M), the bound-state binding energy and 
M the charmonium mass (in order to be consistent, we 
keep only term linear in b, since the 0(q 2 /m 2 ) terms have 
been neglected in the propagator). 



k h /j, 





(b) cc -> a 



FIG. 1: Effective coupling between a cc pair and two photons 
(a) and a dilepton pair (b). 



spin symmetry (HQSS) and to show that the results are 
identical to those of non-relativistic calculations. 

Defining (OIct^cI?/;) = fyM^s^, we have the following 
expression for the amplitude for ip — > it: 



o(47ra em ) 



M, 



4* 



(4) 



from which we obtain the width (neglecting the lepton 
masses) : 



i 



64ir 2 M,, 



dn\M\ 2 = 



3M- 



(5) 



Using M^f 2 = 12 |V> (0)| 2 H|, we recover the non- 
relativistic result of Kwong et al. . The experimental 
value for the leptonic width of the J/ip (r e + e - (J/tp) = 
5.40 ± 0.15 ± 0.07 keV 0) and its mass (3.097 GeV) fixes 
-omitting NLO corrections for now- fj/$ at 410 MeV. 
For the ip', we correspondingly get /,/,' at 279 MeV for 
r e+e - {ip 1 ) = 2.10 ± 0.12 keV and M r = 3.686 GeV. 

Similarly, with (0|c7 M 7 5 c|7/ c ) = «/ 7)c -P M , the amplitude 
for rj c — > 77 is readily obtained: 

u 



M 17 = -AiQ 2 {kKat em ) 



M 2 

Vc 



(6) 



from which we obtain the rj c (lS) width (with b rjc ~ 0) 
1 1 f , o 4nQ 4 c a 2 em f 2 



r 77 (r? c ) 



2 64tt 2 M„ 



dn\MY 



M, 



the factor ^ being the Bose-symmetry factor. 

As suggested by HQSS, let us now suppose the equality 
between fji^ and /^ c , enabling us the following evalua- 
tion, r 77 (?7 c ) = 7.46 keV. 

When NLO corrections are taken into account 1131 . 
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(8) 



with a s — 0.26, r 77 (?y c ) is shifted to 9.66 keV. The latter 
agrees with the world average value 7.4±0.9±2.1 keV Q 
in view of the large statistic and systematic uncertainties 
in the measured value. This indicates that our effective 
Lagrangian approach can also successfully predict the r\ c 
two-photon width. The agreement with experiment also 
suggests that there is no large spin-symmetry breaking 
term in the charm vector and axial- vector current matrix 
elements. We now use the same effective Lagrangian and 
HQSS to compute the r\' c two-photon width. 



III. HEAVY-QUARK SPIN-SYMMETRY 
PREDICTION FOR r 77 (ry c ) 

First, we want to redo the calculation of T^(ip) and 



IV. HQSS PREDICTIONS FOR T 77 (^) 

We now turn to the excited states. Extrapolating 
HQSS to 2S states, i.e. f^i = f n > , and neglecting bind- 

f 



r 77 (r/ c ) through the simple application of heavy-quark in S ener gy effects, we obtain T^tQ 
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3.45 keV, which is more than twice larger than the eval- 
uation by CLEO (1.3±0.6 keV) although nearly in agree- 
ment with Ackleh et al. (3.7 keV), Kim et al. 
(4.44 ± 0.48 keV), Ahmady et al. @ (5.7 ±0.5 ±0.6 keV). 

Binding-energy effects are easily taken into account 
by introducing a correcting factor such that r 77 (7^.) can 
be written as as a function of r 77 (?7 c ), T e + e -(J/ip) and 



e+e" 



- (ip') as follows 



^■yfiVc) ~ L 7 




K M vc 



M^+b^M^I Ml 



Ml 



[J/^Mj/tp, 



This gives 



r 77 (^)=4.1keV, 



(9) 



(10) 



It has to be noted however that the experimental val- 
ues of F 77 (?y c ) are affected by a large systematic uncer- 
tainty related to the branching B(r\ c — > KKir) and the 
evaluation of F 77 (r^,) done by CLEO was realised by as- 
suming B(r) c — v KKt:) = B(rf c — > KKir) which is only 
to hold approximately. This assumption also allows an 
extraction of B(B — » Kr]' c ) from the Belle measurement 
of the ratio (B(B -> Kr)' c )x B(r}' c -> KKtt))/{B{B -> 
i^7? c ) x 6(?7 C -> X^tt) H|. The value of the ratio 
B(B — > Krj' c )/B(B — > if?7 c ) thus obtained seems to agree 
with a theoretical prediction using QCD factorisation 
model for colour-suppressed -B decays with a charmonium 
in the final state 23] . Thus the assumption of the approx- 
imate equality between the rj' c — > KKt: and ?7 C — ► KKir 
branching ratio seems to be justified to some extent. In 
this case, the CLEO low value for the ratio (B(r)' c — > 
77)^(77^ -> KKtt))/(B(t] c -> 77)K(r/ c -> MTtt)) would 
imply the small — > 77 decay rate quoted above. 



therefore the introduction of differences in the mass of 
r] c and rf c increases the discrepancy with the experimen- 
tal result obtained by CLEO. Note that, up to correc- 
tions due to differences in the scale of a s , the radia- 
tive corrections are canceled in Eq. as well as in the 
formula giving the first quoted value, 3.45 keV. If one 
wanted to introduce relativistic corrections in the spirit 
of NRQCD, one would expect them to cancel also, follow- 
ing Eqs (4.3c), (4.3d), (A31), (A32c), (A34) and (A35a) 
of Bodwin et al. El. 



There exist however models that are able to reproduce 
correctly L 77 (^) but, in general, they tend to underesti- 
mate L 77 (^ c ). Indeed, Miinz [lj| predicts 3.5 ± 0.4 keV 
for r] c and 1.38 ± 0.3 keV for r]' c , Chao et al. E3 5.5 keV 
for ij c and 2.1 keV for rf c , and Ebert et al. [l2( 5.5 keV 
for ij c and 1.8 keV for rj' c (see also the results of [25^. 
This clearly points at a specificity not yet understood 
of the rf c decay. All the theoretical predictions and the 
experimental measurements can be found in Tabic (P). 





Experiments 


This paper 


Ackleh [4] 


Kim [5] 


Ahmady [6J 


Miinz fill 


Chao [10] 


Ebert [12] 


Tic 


7.4 ± 0.9 ± 2.1 (PDG [7]) 


7.5 - 10 


4.8 


7.14 ±0.95 


11.8 ±0.8 ±0.6 


3.5 ±0.4 


5.5 


5.5 




1.3 ±0.6 (CLEO [3]) 


3.5-4.5 


3.7 


4.44 ± 0.48 


5.7 ±0.5 ±0.6 


1.38 ±0.3 


2.1 


1.8 



TABLE I: Summary of experimental measurements and theoretical predictions for r 77 (^ c ) and r 77 (^c). (All values are in units 
ofkeV). 



V. CONCLUSION 

Whereas heavy quarkonia are supposed to be reason- 
ably described by non-relativistic approximations, some 
works (e.g. [2(| |23, H3, ) have pointed out that non- 
static effects within quarkonium (especially in radially- 
excited states) should not be neglected without further 
considerations. 

In the case of rj c , we have seen that the simple applica- 
tion of HQSS gives a reasonable estimate of the width 
compared to the world-average experimental measure- 
ments. For rj' cl we have obtained the same discrepancy 
as other models. On the other hand, we have shown here 
that the introduction of binding energy in this calcula- 
tion introduces a correction of about 20 % but worsens 



the comparison with the CLEO measurement. 

When one considers the ratio of the two decay widths, 
radiative corrections cancel out, up to effects due to 
changes in the renormalisation scale. This might slightly 
affect the results, but not sufficiently to recover agree- 
ment with data. Of course, heavy-quark spin symmetry 
(or, equivalently, the equality between the decay con- 
stant for the 3 Si and the 1 Sq) could be broken for excited 
states, but it is quite unlikely that it could be so badly 
broken to explain such a discrepancy. 

Since many other works have shown difficulties to re- 
produce both rjc and rj' c two-photon widths, we are look- 
ing forward for a confirmation of the CLEO measure- 
ments, especially through a better understanding of their 
branching ratio in KKn. 
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